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Abstract: 

We study classical dynamics of an open string tachyon T of unstable Dp-brane cou- 
pled to the gauge field A^. In the vacuum with vanishing potential, V = 0, two 
fluid-like degrees of freedom, string fluid and tachyon matter, survive the tachyon 
condensation. We offer general analysis of the associated Hamiltonian dynamics in 
arbitrary background. The canonical field equations are organized into two sets, fluid 
equations of motion augmented by an integrability condition. We show that a large 
class of motionless and degenerate family of classical solutions exist and represent 
arbitrary transverse distribution of tachyon matter and flux lines. We further test 
their stability by perturbing the fluid equation up to the second order. 



Second half of this note considers possibility of V ^ in the dynamics. We incor- 
porate V in the Hamiltonian equation of motion and consider interaction between 
domain walls and string fluid. During initial phase of tachyon condensation, topo- 
logical defect at T = is shown to attract nearby and parallel flux lines. The final 
state is fundamental strings absorbed and spread in some singular D(p — 1) brane 
soliton. When string fluid is transverse to the domain wall, the latter is known to 
turn into a smooth solution. We point out that a minimal solution of this sort exists 
and saturates a BPS energy bound of fundamental string ending on a D(p — 1) brane. 
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1. String Fluid and Tachyon 



Decay of unstable D-branes [JTJ] have served as useful laboratories of understanding 
off-shell structure of string theory. After several years of intense study, we seem 
to have fairly good idea how lower dimensional, stable D-branes are formed out of 
decaying unstable D-branes. They arises as topological defects, familiar in ordinary 
field theories, although solutions tends to be a bit singular than usual @, || |). Among 
very well understood are how the relevant space-time Ramond-Ramond charges are 
generated via winding numbers and topology of gauge bundles . 

One of reasons that we should anticipate lower dimensional D-branes to emerge 
at the end of day is the charge conservation. Stable D-branes carry Ramond-Ramond 
gauge charge |6[], which can be absorbed in unstable D-brane and transmutes into 
world- volume gauge field configurations Q. Charge conservation in the space-time 
begs for a mechanism for recovering these conserved quantum numbers, and the only 
objects that can carry such charge after the unstable brane annihilate itself, would 
be the stable D-branes themselves. Thus, we must somehow be able to reproduce 
the D-branes within the dynamics of unstable D-branes. 

One would expect the same logic should apply to fundamental strings 0, |J , say, 
to infinitely long semiclassical ones at least. Yet, emergence of fundamental strings 
has proved a much more challenging problem. Some tantalizing hints have been 
accumulating via study of low energy effective action of D-brane decay, nevertheless, 
and in this note, we will give a comprehensive review of the classical low energy 
dynamics with a commonly employed effective action, conceived by many authors, 
and study various aspects with a view toward formation of fundamental string. The 
form of the Lagrangian we will study is of the form [11, 12], 13, [14|, 15 , plj , 

-V(T)y/-Det(ri + F + VX J VX J + VTVT), (1.1) 

where the tachyon T appears inside the determinant on equal footing as the trans- 
verse scalars X 1 . V(T) is everywhere nonnegative and has a runaway behavior 



[[17], |TI| |19j. That is V(T) vanishes exponentially at T = ±oo. For first half of this 
note, we will consider V = strictly. Latter part of the note will incorporate V ^ 
and compute its effect near domain wall formation. 

The earliest evidence that this low energy field theory is quite unconventional 
came in Ref. ||, which considered pure gauge dynamics in decay of unstable D2- 
brane. They found that the electric flux lines become free in that no transverse 
pressure is present and also that the tension density of the flux lines obeys a BPS- 
like property which would be normally seen in fundamental string. This degenerate 
behavior of flux lines were later found to persist in higher dimensional case as well 
2(J, once V(T) is taken to vanish. 

This pressureless collection of flux lines were dubbed "string fluid," in obvious 
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reference to the fact that they carry fundamental string charges.* One outstanding 
question is whether and how properly quantized fundamental string emerges from 
this string fluid. Tantalizingly similarities already exist between string fluid and 
classical fundamental string raising hope that via some confinement mechanism we 
may recover fundamental string from decay of unstable D-branes. Classical properties 
of string fluid are studied in detail in Ref. [^tj . 

It turns out such fluid-like behavior is not limited to the gauge sector but extends 
to the tachyon. When the final state is static and homogeneous, an on-shell condition 
reads gfj |§ [25| 

l = E i E i + f 2 , (1.2) 

where E{ = F 0i are components of the field strength. Energy density of any motion- 
less state is composed of two components 



Hmotionless = \J TT l TT l + 7Tj., (1.3) 

where n l is the conserved electric flux, such that Vj7r* = 0, while 7r T is the canonical 
momenta of T. The kinematics are such that we have a relationship, 

Vttt = Ei/f. (1.4) 

Pure string fluid emerges when T = = 7r^, while the other limit E — involves 
energy density composed solely of ttt- In the latter limit, all energy is carried by 



dust-like matter, known as tachyon matter fl26| , yH 

One of less understood aspects of the combined system of string fluid and the 
tachyon matter is how the two components interact with each other P5| , ^| 

This system was addressed comprehensively and without any approximation, in 
Ref. [25fl. Among explained are the origin of the pressureless nature of these fluid, 



and also exact Hamiltonian, canonical field equations, and energy-momentum tensor 
are written down. The combined fluid shows character of 1+1 dimensional system 
with variable "speed of light", which depends on the composition of the two fluid. 
Although exact "static" solutions were found, generic aspect of the dynamics has been 
poorly understood. In this note, we will address this issue and try to understand 
how one fluid component react to the presence of the other. 

We provide a comprehensive analysis of the classical dynamics in the Hamilto- 
nian formulation, taking into account possible coupling to background metric. One 
purpose is to clarify the relation between the fluid equation of motions and the Hamil- 
tonian equation of motion. Previously, the former has been derived from combination 
of the latter and the energy-momentum conservation. In this note, we will describe 
in what sense solution to the fluid equation gives solutions to the field equations. 
Sections 2,3, and 4 are devoted to this. 

*This fluid- like behavior had been also noticed as the strong coupling limit of Born-Infeld system 

fl2ll I2II2I. 
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In particular, this allows us a rather sweeping characterization of all static (mean- 
ing that physical momenta vanishes identically) classical solutions and shows their 
rather huge degeneracy. After isolating such static solutions, we test their stability 
by perturbing the fluid system in section 5. Perturbation up to 2nd order is per- 
formed, which looks pretty involved, and we argue that all physical effect from such 
perturbation may be understood as a consequence of continuity of the fluid and lacks 
any destabilizing interaction. 

In fact, the aggregation of energy density and flux lines turns out to be tied 
to formation of domain walls, instead. In section 6, we extend the Hamiltonian 
dynamics formulation to include possibility V 7^ 0. We show that the domain- wall 
tends to attract nearby and parallel flux lines via a short range attraction. The 
effective range of this attraction is given by region with finite V, and collapses as 
stable D(p — 1) branes form. However, once the domain wall formation is complete, 
the range of attractive interaction become arbitrarily small, and such aggregation of 
(parallel) string fluid is no longer favored. ^ 

In section 7, we consider another configuration involving string fluid and domain 
wall, where the former lies orthogonal to the latter. Smooth solutions of this kind 
were recently written down, and here we reproduce them from our Hamiltonian 
description above. In particular, we show that a minimal solution exists and saturate 
BPS bound which is normally associated with 1/4 BPS configuration of fundamental 
string(s) ending on or passing through a D-brane. We close with a summary. 



2. Fluid Equation, Integrability Condition, and Classical So- 
lutions 



Fluid equation for V = was first written down in Ref. [p0|| . The main purpose 



there was to understand gauge dynamics, but the analysis did deal with possibility 
of turning on transverse scalars and also dynamical T provided that the kinetic term 
of the latter shows up on equal footing as transverse scalars X 1 . This is precisely 
the Lagrangian in question. Furthermore, any scalar whose kinetic term appears 
this way can be treated as if it is a component of gauge field along some hidden 
direction, and manipulations for gauge field carries over almost verbatim. Because 
of this, the result carries over immediately to the system of string fluid coupled to 
tachyon matter. In this section, we will rewrite the result with tachyon field explicitly 
expressed, and consider its implications. In section 3 and 4, we will elevate this to 
the general background with curved metric. 



^We must caution readers that this effect is unrelated to usual attraction between D-branes 
and parallel fundamental string. The latter arises from exchange of closed strings between the two 
objects, and thus is normally associated with open string one-loop. The net effect is the same, 
nevertheless, so such a tendency to form bound state of a ~D(p — 1) and fundamental strings is not 
too surprising. 
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The main idea is to deal with the dynamics from the Hamiltonian viewpoint. 
Let us introduce extended gauge field Am = (A^At) so that At = T. Then the 



Lagrangian of tachyon coupled to a gauge field may be written succinctly as [30 



-V{T) y/-det{riMN + F MN ), (2.1) 

where F MN is the field strength associate with A M with the understanding that 
dx = 0. Since we will be using Hamiltonian formulation we will be separating 
out time direction from the "spatial" ones. For the latter we will use subscript 
m = l,---,p,T and i = 1, •■•,£>. Using Legendre transformation, we obtain the 
following Hamiltonian of this system 

H = H- A d t ir\ (2.2) 

where 

H = y/ir m n m + P m P m + V 2 det(5 mn + F mn ). (2.3) 



van- 



Here 7r m is the canonical conjugate momenta, defined as dC/ dA m , while P m = 7r n F„ 
Pi with i = l,2,---,p correspond to the conserved momenta associated with the 
translational symmetry. 

In section 2 through section 5, we will consider the limit, V = 0, which repre- 
sent the final stage of tachyon condensation where tachyon T is rolling away to oo 
everywhere. The Hamiltonian (in temporal gauge) is exceedingly simple, and has 



H = y/7T m 7C m + P m P m . (2.4) 

Note that despite the vanishing Lagrangian as V — > 0, the Hamiltonian remains 
finite. 

2.1 Fluid Equations 

Half of the canonical equation of motion 

5 



5A r 



dx p H, (2.5) 



combined with conservation of energy momentum gives immediately the following 
fluid equations |pOfl 

d n m + v'din 171 = n'div" 1 , 

d v m + v%v m = n'din" 1 , (2.6) 
where the vector fields n and v are defined as 

7r = Hn, 

P = Hv, (2.7) 
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which satisfies the constraints 




n m v m = 0. 



(2.8) 



The evolution of the energy density H is then determined via, 



d H + di (Hv l ) = 0, 



(2.9) 



and finally 7r l must satisfy Gauss's constraint, 



diir 1 = 0. 



(2.10) 



The Noether momenta Pi has the following simple expression 



Pi = -Fijit 3 - diT7r T = Tr m F, 



(2.11) 



The last component of the vector P, P?, is not really a conserved momenta since T-th 
direction is a mere mathematical device of convenience. It is nevertheless computed 
by pretending T-th direction exists as a translationally invariant spatial direction, 
and thus given by the combination 



Most of quantities here have simple physical interpretation, is nothing but con- 
served electric flux, while nr is the tachyon matter density. H and Pi are conserved 
energy and momentum density, respectively. The last quantity Pt measure inho- 
mogeneity of tachyon T along the flux direction. The fact that this appears in the 
Hamiltonian separately implies a rather anisotropic behavior the system. 

2.2 Integrability Conditions 

Although fluid equations are self-contained, a solution may not give automatically a 
solution to the original field equations. This is because that the fluid variables above 
are naturally formed from canonical variables and does not produce elementary fields 
A m directly. A further set of first-order equations must be solved. These arise from 



and are exactly half of the Hamiltonian equations of motion. 

These equations are intimately related to the fact that the determinant part of 
the Lagrangian vanishes. The Hamiltonian with V kept is such that 



P T = d{Tn l = n m F mT . 



(2.12) 




(2.13) 




(2.14) 



for some matrix C independent of ir m . Then the Lagrangian is 



L = A m n m -H = n m ( [ h) - H = 0(V 2 ). 




(2.15) 
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On the other hand the Lagrangian is of the form 

-Vy/-det(---), (2.16) 
so the on-shell value of Lagrangian must be such that 



-det(- • •) ~ V -»• 0, (2.17) 

which is precisely the condition that generalizes, 

1 - T 2 - EtEi -> 0, (2.18) 

of the homogeneously rolling tachyon. 

In any case, these equations can be expressed as an algebraic equation for the 
field strength associated with elementary fields A q in terms of fluid variables, H, n, 

v; 

E m f^m Pmn^ni (2-19) 

where E m = F 0p . This, combined with the Bianchi identity 

Fran d ni E n d n E m , (2.20) 

gives an evolution equation for F mn . 

On the surface, thus, it may seem that we have split the system into two steps; 
solve the first-order fluid equation and then solve another first order equations for 
elementary fields in the background of H, n, and v. However, there is one subtlety 
here. The two fluid variables v and n are algebraically related as 

^Vn. UnF-rimi (2.21) 

and there is a logical possibility that one may not find such F mn as a solution. Only 
if there is a solution F mn consistent with ( |2.21 ), the solution to fluid equation would 
be acceptable. 

In the simplest case of pure tachyon, this integrability condition appears in a 
particularly simple manner. The only elementary field is At = T in that case, and 
we have = 0. Then (|2.21|) states that 

Vi = -diTn T = T , =, (2.22) 

while ( |2.19| ) gives, 



f = (1 + (djT) 2 ) n T = ±y/l + {d t Ty. (2.23) 

In this case the integrability condition may be imposed as a familiar local condition 
on the fluid system as 

*te)- % te)=°- (2 ' 24) 

In section 3, we will further consider the pure tachyon case and reformulate these 
conditions in a manifestly relativistic manner. 
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2.3 "Static" Solutions 

Let us characterize all static solutions. By "static" we mean absence of physical 
momentum, Pj = 0, and this immediately constrains us to look for configurations 
with Vi = 0. One of the algebraic constraint then implies that ti t vt = 0. This can 
be achieved by v t = in addition to v j = 0, and if we choose this,* fluid equation 
collapse to 

d H = 0, (2.25) 
d n m = 0, (2.26) 

and 

n%n m = o, (2.27) 
di (Hrt) = 0. (2.28) 

Because of ( [2.271 ), A ux lines should be straight, and we may as well associate its 
direction with x 1 while denoting the remaining p — 1 directions by x 2,3 '"' ,p . The most 
general "static" solutions (with vt = 0) to the fluid equations are time-independent 
and ^-independent distribution of flux n 1 and tachyon matter tt t . By construction 
the same property holds for H. Then, all "static" classical solutions are characterized 
by the following two arbitrary and independent (nonnegative) functions^ 

together with the choice of x 1 direction or equivalently the choice of the electric flux 
direction. The energy density is 

H(x 2 , ...,x p ) = ^ (tt 1 ^ 2 , x 3 , • • • , xp)) 2 + (tt t (x 2 , x 3 , ■ ■ ■ , xp)) 2 . (2.29) 

In section 5, we will analyze classical stability of this large family of static solutions. 

We still need to test whether this "static" solution satisfies all the integrability 
condition. Namely, we need to find a gauge field, A m , solving (|2.19|) , such that ( |2.21| ) 
holds. Since v = by construction, we have 

E m = n m = 7i m /H, (2.30) 

J The other choice of tit — leads to a different class of solution that involves domain walls 
pi] , |32| , |33| . Since 7^ v t — fiidiT implies nonvanishing gradient of T along the flux direction, 
such a configuration will gave T = somewhere along the flux line direction. The final state would 
corresponds to a domain wall threaded by a transverse string fluid. Because of this, one can no 
longer restrict to V — > limit, and one must solve for the full equation of motion. See section 7 for 
this class of solutions. 

§An exceptional case is when it 1 also happens to vanish; tt t could then be an arbitrary function 
of 3.1,2,3,- -, P _ 
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which has a solution in the temporal gauge Aq = 

A m = tn m (x 2 ,---,x p ). (2.31) 

Then the magnetic field is 

-Ftj =Fji = tdjUx, 

-F T j = F jT = tdjflT, 



vanishing v via v m = n n F nm} we find 



for j — 2, 3, • • • , p. All other components vanish. Testing whether this gives back a 
«i = 0, 

^ = -Fj-i/n - F iT n T = -tdj I J = 0, 

w T = 0. (2.32) 

Thus our "static" solution is integrable and thus acceptable as a solution to the 
original field equations. 

3. Tachyon Matter in General Background 

The tachyon effective action on an unstable Dp brane in general gravity background 
is given by 

S = - J d p+1 xV (T) y/^g s/ 1 + g^dpTdvT, (3.1) 

where g = det g^ v and V(T) is the tachyon potential, and we defined the general 
metric as 

ds 2 = g^dx^dx" = -N 2 dt 2 + h ij (dx i + L i dt)(dx j + Vdt) 

= sTd^d v = - m) 2 + Wdidj, (l < i,j < p). (3.2) 

This general metric is determined by the lapse function N, shift vector L l , and spatial 
metric hij of the p-dimensional hyper-surface at a given time slice. Using this general 
metric (|3.2j), we can rewrite the action (|3.1| ) as 



S = J dt J d p xC, 

C = -V(T)NVhVx~, (3.3) 

with h = det h^, g = —N 2 h and 

X = 1 + grdpTdvT 

= l-j^(f- Vd{T) 2 + WdiTdjT, (3.4) 



- 9 - 



where the doted notation denotes time derivative. Then the Hamiltonian is obtained 
by 

H = vrT - C 

= 7if + V{T)Ny/h\fk 



= N y ir 2 + (irdiT^iwdiT) + V 2 h(l + WdJTdjT) + nVdiT, (3.5) 
in which we have defined the conjugate momentum as 



9T NVX 



(3.6) 



After tachyon condensation, i.e., in V — » limit, the Hamiltonian equations are given 
by 



f = — = NJl + hvdiTdjT + UdiT, 
On 



&H 



TV 



d, 



Nirh ij djT 



dT ~'° \ ^l + h kl d k Td{T , 
Let us consider the following Lorentz invariant matter density 



(3.7) 
(3.8) 



p a d a T 



where we define 



V 

/x 

V 

s/X 



N 2 



(T - L%T) 2 - h ij diTdjT 



dC 



d(d a ry 

so that p° = n T . In the last step of the above equation, we used the Eq. 
Using the equations (|3.7p , ( p.8|) , we obtain 



(3.9) 
(3.10) 



-g 



-g 
l 



M-g 



'-g 



-A 



\N 2 



-d 7i + d t {nU) + $ 



0. (3.11) 



We used the Eq. ( |3.8| ) in the last step of the Eq. ( |3.11| ). Then the energy momentum 
conservation implies that 



= V M T, 



V 



=duTd v T 



fid^TWd u T. 



(3.12) 
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All of these have a rather obvious interpretation once we identify fi as the invariant 
matter density and — <9 M T as the velocity field PH . Equations (|3.11| ) and ( |3.12j ) 
then implies, 

V Q (/if/ a ) = 0, (3.13) 
U a V a U» = 0, (3.14) 

which are nothing but the continuity equation and the geodesic equation for dust. 
These are the two fluid equation for case of pure tachyon. 

If we were treating the density fi and velocity field as elementary quantities, 



the fluid equations (|3.14|) must be augmented by an integrability condition 

- VuU^ = 0, (3.15) 

to make contact with the original field theory. This of course gives = —d^T for 
some function T, and we would interpret this T as the original tachyon of the system. 
Also recoverable from the field equations is the fact that 

= -1, (3.16) 

which is a kinematical constraint, saying that the velocity field = —<9 M T arises 
from affine parameterization of the geodesies. Obviously it does not affect motion of 
the tachyon matter, and simply gives how T should be solved for, given any particular 
set of trajectories of tachyon matter. 

Thus the tachyon matter in the tachyon condensation limit corresponds with an 
ideal (rotationless) fluid, moving freely along geodesies, with the trajectory being 
afnne-parameterized. In particular, this implies that the tachyon matter clusters 
under the gravitational interactions just as ordinary matter does. 

4. String Fluid Coupled to Tachyon Matter 

Here we repeat the above analysis with gauge field included. The fluid equation for 
string fluid were first written in Ref . |2(| , and here we show their completeness again 



in that all nontrivial dynamical information can be recovered just from solving the 
fluid equations. The effective action for tachyon coupled to Abelian gauge field on 
Dp-brane can be written by 



with 



S = J d p+1 xC, (4.1) 
C = -V(T)V^X 

X = det{g^ + d^Td v T + F^), (jm, v = 0, 1, • • • ,p), (4.2) 
F = d^A v - d v A^ 
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where g^ is defined in Eq. ( j3.2|) . 

In calculating the determinant in Eq. Q4.2|) , it is convenient to consider the 
tachyon field T as (p + l)-component of the gauge field, as we already used in section 
2. The metric is also extended to include a flat fictitious direction x T . Then the 
Eq. ( |4.2| ) can be rewritten by 

X = det(g MN + F MN ), (M,N = 0,l,---,p,T), (4.3) 

where F^t — d^T since dx = 0. Using Legendre transformation, we obtain the 
following Hamiltonian of this system 

H = Ny/n m h mn n n + P m h mn P n + V 2 det(X mn ) - P m L m - Aody, (4.4) 

where n m = SC/SA m , P m = 7r n F nm . Canonical field equations for this system can be 
found in Ref. p5f , and in the following we will reexpress these in terms of fluid-like 
variables and also generalize it to include nontrivial background metric. 

4.1 Generalities: Case of an Unstable D2-Brane 

Our purpose in this paper is to understand the interaction between tachyon matter 
and electric flux line. To accomplish this purpose, we concentrate on the simplest 
case, D2-brane. We believe that the tachyon condensation on D2-brane contains all 
nontrivial characteristics of the interaction of tachyon and flux in general Dp-brane 
decay. The determinant in Eq. (f4.2|) on D2-brane system is given by 

X = g ( 1 + V - -j-(F*F) 2 ) = gX, (4.5) 



2 16 



with 



F 2 = F MN F MN , F*F = F* MN F MN , F* MN = -±— e MNPQ F PQ , (4.6) 

2 v 9 

where M, N = 0, 1, 2, T and ^ MNp Q is the Levi-Civita symbol and we choose eoi2T = 
1, e 012T = — 1. Energy-momentum tensor T^ v is^ 

T „ s * « = (4.7) 
yf^gog^ \JX 9 

where C MN is the symmetric part of the cofactor Cmn for Xmn = 9mn + Fmn which 
is given by 

C MN = g fMN^ + 1 F 2 } + p MN + F MPp p N _ ^ F *MN^ p) \ (4g) 



^Bear in mind that, since we artificially introduced x T direction, there is no conserved momentum 
associated with this direction. 
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Equation of motion for the gauge field Am and conservation of energy-momentum 
are expressed by 




(4.9) 
(4.10) 



where C^ N is the anti-symmetric part of the cofactor in Eq. ( |4.8|) . Conjugate mo- 
mentum for gauge field and conserved Noether momentum which is induced by in- 
variance of the spatial translation are given by 



7T 



V C / 1 



P 



T/ r<0m 

A TP 



(4.11) 



In the next sections, we will do that in tachyon condensation limit in flat and curved 
space. 

4.2 Fluid Equations in Flat Space 

In flat space metric tjmn = diag(— 1, 1, 1, 1), the Hamiltonian density ( |2.2|) in tachyon 
condensation limit (V = 0) is given by 



H = \Jix 2 + (B x vf) 2 = -^=(1 + B 2 ), 

vx 

where we use A = gauge and define E m = F 0m , B m = e mn iF nl and thus 

X = 1 + B 2 - E 2 - (E ■ Bf = 0, 



with an explicit expression 

E = (Ei, E 2 , T), 

With this, we find 



B = (d 2 T,-d 1 T,B). 



(4.12) 



(4.13) 



(4.14) 



V 



Tin 



E m + B m (E ■ B) 



[B X 7?) 



V 



B n E h (4.15) 



and under the condition (|4.13| ), 

V 



^jrnn 



'X 
V 



: °5 



P jr — nr P 

H ' 

PjPj - TTjTTj 
H ' 



(4.16) 
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Equations of motion gives 



d ■ {Hn) = 0, (4.17) 
ft + (v ■ d)n = (n • d)v, (4-18) 

when augmented with the energy conservation, 

H + d-(Hv)=0. (4.19) 

The momentum conservation may be written as 

{/+ (v-d)v = (n-d)n } (4.20) 

where energy conservation is taken into account. Actually, the last component of 
( [4. 20| ) does not arise from the conservation law, since no such conservation law exists 
for Pt- Rather it is a linearly dependent equation that may be derived from the rest. 
However, we include this time evolution equation for vt for notational convenience. 
As already described in section 2, the vectors n and v are defined as 

7T P 

n=-, v = -. (4.21) 

The two vectors, n and v satisfy the constraints 

n 2 + v 2 = l, n-v = 0. (4.22) 

Thus, the dynamics produces a set of self-contained first-order fluid equations for H 
and n and v. Although these equations are more natural in Hamiltonian formula- 
tion, we stick to Lagrangian formulation because the latter is more susceptible to 
incorporation of curved background. 

4.3 Fluid Equations in Curved Space 

In the tachyon condensation limit (V = 0), the Hamiltonian in curved space is written 
by 

H = N^Ti m h mn TT n + P m h™P n - L m P m , (4.23) 

where we use A = gauge and the X is defined in Eq. ( fO| ) and to ensure the 
finiteness of Hamiltonian density we have to set X = 0. After some calculation, we 
obtain the following two relations which is similar to the flat space case 

V pm n _ m pn 
r /^min _ 



/ X 

V pmpn _ m n 

—Cs n = = Wo , (4-24) 
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where T 00 is the (OO)-component of the energy-momentum tensor defined in Eq. (|4.7| ) 
and we define the upper indexed P l as 



P m = g m0 7T n F n0 + g mn P n . (4.25) 



Using the relation in Eq. ( [4.24| ), we can rewrite the equation of motion which is 



expressed in Eqs. (f4.9|) as follows : 



7T 



Vj ( —^= ) = 0, (y = case), 



where 



V ( j + Vj f ^_ ^ ) = 0, (z/ = *case), (4.26) 



H = ^T m = ^-j= v^W^ + P m h mn P n . (4.27) 



And the conservation of energy-momentum (|4.10|) is given by 



f H \ P i 
V ( -5= ) + V f -= = (z/ = 0case), 

pi \ /pi pj _ i j \ 
V ( -= + Vj ^_ rr ) = (u = i case). (4.28) 



This is an analog of ( |4 . 1 7| - ^720D in curved spacetime. 

The form of fluid equations take particularly simple form when the shift vectors 
happen to vanish, IS = and N is constant. This would be the case when we are 
considering a cosmological scenario, for instance. With this, Vo effectively collapses 
to (1/N)d t , and the Eqs. ( |4.26| ), ( |4.28| ) are summarized by 



V' • (Hn) = 0, (4.29) 

H + V' • (Hv) = 0, (4.30) 

n+{v-V)n = {n-V')v, (4.31) 

v + (v- V')v = (n- V')n, (4.32) 

where V = (V^, V 2 , 0) is the covariant derivative defined for the spatial part of the 
metric, h mn , and the dot represent the time derivative (1/N)d t . We have also defined 

n m ^ ~v m = (4.33) 

with constraints 

h 2 + v 2 = l, n-v = 0. (4.34) 
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4.4 Integrability 

Since the quantities that enter the definition of TC, n, and v, are naturally canonical 
variables, solving the first-order fluid equation will not generate the elementary fields 
T and directly. For these, we have to solve for another set of first order differential 
equations. 

We could have derived the above fluid equation from the canonical formulation, 
and if we did so, half of the canonical equations would remain unused and simply 
relate conjugate momenta to time-derivative of elementary fields. At least when we 
make use of energy-momentum conservation directly. In terms of the above variables, 
the remaining equations may be written as 

E = n-B xv. (4.35) 

This determines the evolution of the "magnetic field" B via Bianchi identity as 

B = d x n-d x (B x v). (4.36) 

Thus, B may be solved for after v and ft are determined. When this is done, E is 



also determined via ( 4.35 ). 



However, since B enters the fluid degrees of freedom via the identity, 

v = Bxn, (4.37) 

we need to make sure that this does not further restrict solutions to the fluid equations 
alone. A priori, it is unclear whether for all solutions to the fluid equations we can 
find B such that it solves the evolution equation and is consistent with this algebraic 
constraint as well. Only if such a B exist for the solution to the fluid equation, we 
can say that the solution is physical. 

In the curved background these integrability conditions are expressed as 

E m = ^= (h mn h n + F mn v n ) - F mn L n , (4.38) 

h mn v n = F mn n n . (4.39) 
plus the Bianchi identity 

F = d E - d E (4 40) 



5. Stability Analysis of Classical Solutions 

In this section, we perform stability analysis of "static" solutions of section 2. Be- 
cause of lack of transverse pressure, dimensionality of the unstable D-brane does 
not really matter, so we will confine our analysis to the case of unstable D2-brane. 
The authors of Ref. recently debated against the assertion that the degeneracy 
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implies stability. Here we demonstrate the stability explicitly by solving dynamical 
equations of motion, and thereby invalidate the criticism. 

Two cartesian coordinates are denoted as x and y. We will solve the Eq. ( |4.17| )- 
( |4.20| ) in flat space using perturbation theory. The Eqs. ( ^4.18|) , ( f4.20|) may be con- 
veniently rewritten as 

a-(b-d)a = 0, b+(a-d)b = 0, (5.1) 

where we define 

a = n + v, b = n — v. (5-2) 
The constraints ( |4.22j ) become 

a 2 = 1, b 2 = 1. (5.3) 

Let us then consider the following small fluctuations around the background fields 
ao, bo which is given by the "static" classical solutions given in section 2. This means 
that Sq = bo = uq with 

n = (n<te(y),0,7ior(y))- (5.4) 
Keep in mind here that < (n T) 2 = 1 — (^ox) 2 < 1- 

5.1 First Order 

Expanding around such a solution 

a = a + cr- 1 ' + a*- 2 ** , • • • , 

b = bo + b {l) +b {2) ,---, (5.5) 
and using the Eqs. (|5.1| ), fl5.5|) , we obtain the first order perturbation equations 

2n 0a; y 

ft T?(°) 

5 + 6« = -^-a« (5.6) 

where we define 

x ± = x ± n 0x t, d± = ±-—(d t ± n 0x d x ). 

2n 0x 

The y-components of the first order perturbation equations in Eq. (|5.6|) are expressed 
by two homogeneous first order differential equations 

&_o« = 0, «9+6« = 0, (5.7) 

and the solutions are 



a 



(i) 



f(x+,y), b®=f{x-,y), (5.8) 
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where /(/) is an arbitrary function with arguments x + (x~~) and y. 
Substituting the solutions for ay and into the x-components of the Eq. 
we get 

d-a^ = -M^/V.y), fll^W = -^/(* + ,y), (5.9) 
2n 0a: 2n x 

and the solutions are 



a« = 0(3+, y) - ^ifyr, y), 6« = ftaT, y) - ^F(x + , y), (5.10) 
where g and g are arbitrary functions, = d y n 0x , and we define 

/X + fX~ 
dw f(w,y), F(x~,y) = / dw f(w,y). (5.11) 

Using Eq. (|5.9|) , we obtain free wave equations, d + d~a x : (&£ ) = 0, i.e., 

- n?> = 0, (d 2 - n ^ 2 ) = 0, (5.12) 

this means that there are propagating modes along background flux line |25|, and 
no net effects which deforms the distribution of the background flux line and tachyon 
matter in the first other perturbation. 

5.2 Second Order: Low Frequency Limit 

Now let us consider the second order perturbation of the Eq. ( |5.1| ). General for- 
malisms for the second order perturbation are analyzed in Appendix C. As a mean- 
ingful choice for the solution of the first order perturbation equations to investigate 
the second order ones, we consider the following configurations 

f(x + , y) = u(y) + c(x + , y), (5.13) 
Rx~, y) = ~u(y) + c(x~,y), (5.14) 

g(x + , y) = --^u(y)x + + d(x+, y), (5.15) 
2n 0x 

Tl! ~ 

g{%~, V) = 7T^u{y)x' + d(x~, y), (5.16) 
2n 0x 

where c(c) and d(d) are arbitrary oscillatory functions. In other words, we allow a 
net velocity u(y) along y-direction in the first order perturbation. 

Since the second order perturbation is pretty involved, let us take some simpli- 
fying limit. We could for instance consider taking a very low frequency limit. After 
all, if there is a confining or dispersive effect, it should show up here. In the current 
general setup, we may achieve this by taking the limit 

/SrV.lO-O, ( 5 - 17 ) 
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where f osc represent oscillatory piece and m, n are arbitrary positive integers. We 
use the notation, 

A {m > n \x, y) = =—A(x, y). (5.18) 

In addition, we average over characteristic time-scale of the oscillatory pieces, mean- 
ing that we are mainly interested in net effect rather than exact time evolution of 
the system. We achieve this by averaging the equation over time where in effect we 
set 

1 ft 

(fosc) = Jim - / dwf osc (w,y) -> 0. (5.19) 

L— >oo L I L 

In this limit, the Eqs. ( p.!2|) - ( |0.15|) are reduced to 

(d 2 - n Qx 2 d 2 x )nf = 2uu'n' 0x + u 2 nl, (5.20) 

(d 2 - noM)^ = 0, (5.21) 

(d 2 - n 0x 2 d 2 x )n^ = 0, (5.22) 

(d 2 - noSBfiv® = 0. (5.23) 

Left-hand-sides has the Klein-Gordon kinetic operator for an one-dimensional system 
with "speed of light" equal to a/1 — nl T = \no x \. This is identical to the first order 
case, and implies that small fluctuation would move freely up and down along the 
flux lines, provided that there is no term on the right hand side. 

On average, there is no net force on the second order fluctuation, generated from 



the first order fluctuation, except the right hand side of (|5.20Q . However, this term 
represents a rather trivial effect. u(y) ^ implies that the entire configuration is 
drifting along y direction. Since the flux lines and tachyon matter are distributed 
nontrivially along y direction, this motion will generate time-dependent change of 
n. To see what effect it has, we need to solve for induced by u(y) first. For 
x-independent perturbation, we have 

{7+ (v-d)v = 0, (5.24) 

which gives 

vW(y,t)+vW(y,t) = u(y)-u(y)u'(y)t. (5.25) 

With this velocity field, we may ask how n(y, t) drifts with time. If the only physical 
effect is the drift, the n(y,t) will be identical to n(y, 0) = n(y) where 

y = y+ [ d8v v (f( 8 ), 8 ), (5.26) 

with f(s) is the y-trajectory between /(0) = y and f(t) = y due to the velocity field 
v y (y,t). Despite somewhat involved formulae so far, the relationship between y and 
y is deceptively simple, 

y = y + u(y)t + 0(t 3 ), (5.27) 
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or equivalently 

y = y- u(y)t + u{y)u\y)t 2 + 0(t 3 ). (5.28) 
Thus we find simple drift of the configuration along y-direction gives, 

n(y,t) = n{y) 

= n(y) - t (u(y)d y n(y)) + *- (2u(y)u'(y)d y n(y) + u(y) 2 d 2 n(y)) , (5.29) 



up to 2nd order in time t. Thus the nontrivial term on the right hand side of ( 5.2C|) 
simply represent this drift effect." 

5.3 Second Order: Case of Interlocked Distribution of the Two Fluid 

One large subset of classical solutions we could consider in more detail is those with 
uq x constant. Since n§ x is the ratio between the flux energy density and the total 
energy density, such a solution corresponds to an arbitrary distribution of the energy 
density along y direction, while maintaining the ratio ttq x /ttqx fixed. In this case, all 
terms with derivative on uq x die away, so we have 

d + d^ = --!-fg™ - ^-§ g ™ (5.30) 
d + d-b® = ~fg^ - ^-ggM>, (5.31) 

ZTT-Ox ^flQ x 

d+d-a® = ~ff^ - ^-9f {2fi) , (5.32) 
y 2n 0x 2n 0x 

d + d-b® = L//(i.D _ — 2/( 2 '°). (5.33) 

y 2n 0x 2n 0x 



Since we saw the effect of u(y) is an overall drift of the system along y direction, we 
could safely turn it off in the first order perturbation, /, g, f, and g. Then the first 
two are oscillatory function of x + while the latter two are oscillatory functions of x~ . 
All terms on the right hand side are of the form, 

h(x + ) x h(x~), (5.34) 

for a pair of some oscillatory functions h and h. Such combination of force terms on 
the right hand side cannot generate a net effect, when averaged over time, since the 
resonance effect cannot occur. The right hand sides will drive some oscillation of a 
and b and do not lead to any instability. 

Thus, we conclude that this large class of classical solutions are all stable under 
perturbation and any possible fluctuations move freely along a direction set by electric 



"We could have started perturbation after suppressing oscillatory pieces completely, instead of 
averaging over it, to see such drift effect. If we did that, we would have found exactly ( 5.29| ) as the 
solution. 
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V(T) 




Figure 1: A prototypical form of the potential V(T) as function of T, V(T) 
1/coshT ||, | 



37|, 38, p9fl. We will use this potential for plot of the time 



evolution of flux later in this section. 



flux lines. The origin of this one-dimensional behavior was previously explained 
in terms of collapse of an effective-causal-structure. We again emphasize that the 
integrability condition can at most restrict acceptable solution to the fluid equations, 
so stability under the latter is sufficient to argue the stability under the full field 
equation. 



6. Hamiltonian Dynamics with Potential V 

Much of what we studied above concern a region where tachyon rolls to one side of 
potential. For example we are imagining that T — > oo everywhere. In this section, 
we take into account possibility of V ^ in the dynamics. The formulation here 
should serve useful tool for understanding initial stage of tachyon condensation. 



6.1 Canonical Field Equations with V 

Let us consider the dynamics with V ^ somewhere. For simplicity, we will consider 
unstable D2 brane case again, with all transverse scalars suppressed. Recall that we 
are using the notation introduced in section 2, where a fictitious direction x T is 
employed and T is treated as if it is a component of the gauge field along x T . For 
more details we refers the reader to Eqs. ( 4.14 ), ( 4.1 5|) . 



Half of the Hamiltonian equations of motion 



E = ±-(«-BxP), (6.1) 
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generates the evolution equation for B = (d y T, —d x T, B) when combined with the 
Bianchi identity 

B = d x E. (6.2) 

It is important to note that the energy density H has a V 2 term inside the square 
root, 



H= ^n 2 + P 2 + V 2 (l + B 2 ). (6.3) 

The other half gives evolution equations for electric flux n x y and conjugate momen- 
tum for tachyon 7Tt, 

x, + fl,( '^-'^ + V ' F « )=0, (6.4, 

with V = <9V/«9T and F y = Bey. 

These evolution equations should be consistent with energy-momentum conser- 
vation which now takes the modified form 

H + = 0, (6.6) 
Pl + fl, ( P ^-" V -/ 2 ^ + ^ i )) = 0, (6.7) 

which is almost identical to the previous case except for V 2 terms. We use the 
following facts which are the flat metric version of Eq. 

T oo = V_ c ov = ^ 
VI 

rpOi £jOi pi 

\/x s 



V nik P i P k -n i 7T k - V 2 {5 ik + B i B k ^ 



s 



H 



In all of above equations, it is quite clear that the fluid-like behavior of V — regime 
will be ruined by the potential, and its effect is of order V 2 . 

6.2 Flux Motion near Domain Wall Formation 

In tachyon condensation, V ^ can survive the decay process if there is a topological 
defect. In case of single unstable Dp branes, only possible topological defect would 
be domain walls, separating a region of T = oo from that of T = — oo, and at the end 
of day become D(p — 1) branes in the context of type II theories. We will concentrate 
on initial configuration which will lead to a single flat D(p — 1) brane. With such 
initial configuration, we would like to ask how flux behaves where V ^ 0. 
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Figure 2: Plot of 7f x as function of y, with V(T) = 1/coshT and T = y. 
It shows short range attraction of fluxes toward T = during domain wall 
formation. At late time, T = uy with u — > oo, so the range of the attractive 
force infinitesimal. The final configuration at u = oo corresponds to BPS 
D(p — 1) brane with some fundamental string flux trapped. 



To answer this question, we consider an initial configuration of T such that it has 
spatial variation along y direction and vanishes at y — 0. In addition we assume static 
initial condition, so that neither the tachyon matter nor the string fluid has initial 
velocity. Finally for the sake of simplicity, we further assume a uniform distribution 
of flux lines, lined up along x direction. This can be summarized by the following 
initial conditions, 

T = T(y), 7f = (Vox, 0,0), B = 0, (6.8) 

where tiq x is a constant. 

Now let us consider evolution of Hi, (i — x, y) as follows, 

TTi(t) = TTio + t-ki\t = Q + it 2 7fi| i=0 H . (6.9) 

Under the initial condition ( |6.8j ), the first non-trivial variation of the flux appears in 
TVi\t=o term in Eq. (|6.9j), i.e., 7Tj| i=0 = 0, and is given by** 

**Note that at t = the energy density H reduces to 



while 

B 



v /^ + y2(1 + (W)i 

-MIL' ^= ' ^=^(¥) t=0 - (O0) 
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7T X = — 8y 



7T X Py + V 2 B 

H 



t=o 



t=0 



while Tiy is not generated up to this order, 



7T„ 



0. 



7T„ 



0. 



Thus, we find the leading order time-variation of flux lines is summarized as 



7T, 



7T„ 



1 2 

+ 2* 



H v \H ) H y \H 



+ 0{t s ). 



(6.11) 



(6.12) 



(6.13) 



t=o 



From this expression, it is not too difficult to show that initially uniform n x will 
get redistributed such that flux in the region V ~ 1 tends to gather toward T = 0. 
Since this effect is of order V 2 , which is exponentially small far away from T = 
domain walls, the interaction is of short range. The effective range of this attraction 
is determined by values of V and therefore by the gradient of T along y. With T = uy 
and in unit where a' ~ 1, the range is roughly by ~ 1/u. Thus, the effect is maximal 
when the tachyon begins to roll and die away quickly as the tachyon condensation 
progresses; To reach final stationary state, u has to be infinite, as is well known. 

To illustrate this effect graphically, we draw two figures. Figure 2 is the plot of 
tt x as function of y, assuming an initially linear gradient of T = y. This clearly shows 
local attraction of fluxes toward T = 0. Final figure is a plot of tt x /tc x at y = as 
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a function of ir x /V at y = 0, which shows that the attraction is universal and tends 
to get weakened when flux energy dominates over the tachyon energy. This effect 



explains the phenomenon observed in Ref. JJIJ and produces DQo — 1) branes and 
fundamental strings. 

6.3 Confined Flux Strings are Exceptional and Rare Final State 

Some of more recent proposals for fundamental string formation utilizes configura- 
tions involving T = 0, which would be center of domain wall separating regions of 



T — > oo and of T — > — oo jH], [3o|, [IT]]. In the above, we analyzed the classical dy- 
namics in such situations, and gave a quantitative description on aggregation of flux 
lines near domain walls with T = at the center. We should note here that any such 
process will involve bound states of fundamental string with D(p — 1) branes as the 
final product, and the fundamental string flux will be confined to a co-dimension-one 
hyper-surface rather than a thin, string-like, flux bundle. 

We must emphasize here that the real question lies not in whether there are 
confined flux string but rather how such confined configuration become generic con- 
figuration. In fact, classical solution (even with V = everywhere) that behaves 



exactly like a Nambu-Goto string has been introduced in Ref. |2(J and further stud- 



ied in Ref. pSfl . This classical solution admits infinitely many marginal deformation 
that disperse the flux lines, which makes it very unlikely configuration to form, to 
begin with. 

Decay processes involving topological defects are intriguing in that it could lead 
to this Nambu-Goto string by fine-tuning initial configurations. For instance, one 
could imagine a infinitely long cylinder, R x S p ~ 2 , of a domain wall to which nearby 
electric flux get attached. Depending one details of the initial condition, it is possible 
to imagine that S p ~ 2 part of the cylindrical domain wall collapses after tachyon 
condensation completes and perhaps push all electric flux on its world-volume to an 
infinitesimally thin string. Such configuration will also obey a Nambu-Goto dynamics 
classically |4l|. However, the trouble is simply that, in the phase space of this theory, 
such string-like configurations are very exceptional ones. One must still understand 
exactly how this huge degeneracy is lifted, be it from higher derivative correction or 
via some quantum effect |8|, |44j . 



7. BPS Limit and Strings Orthogonal to ~D(p — 1) Branes 

Fundamental strings orthogonal to a D(p — 1) brane preserves 1/2 of supersymmetry 
left unbroken by the D-brane, so the system of such fundamental strings and a D- 
brane should preserves 1/4 of 32 supersymmetry in spacetime. With the domain 
wall incorporated into the system, we could ask whether there is such a BPS-like 
solution involving string fluid transverse to D(p — 1) branes. Generic solutions with 
right spacetime charges was written down in Ref. PT . 
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This smoothness of the domain wall solution may be contrasted to the singular 
solutions considered in Ref. Q. This smoothness of the solution caution us against 
in using singular domain walls in the presence of transverse string fluid. It was 
recently argued that string fluid would be repelled by D(p — 1) brane, if the former 
try to impinge upon the latter E|. This is based on energetics of the string fluid 



assuming a singular domain wall background. The smoothness of the solution above 
is generically attributable to the presence of transverse string fluid, which shows that 
backreaction due to string fluid may be important in some cases. 

In this section, we will rediscover these smooth solution via our formulation. In 
particular, we find that a minimal solution of a sort exists and saturate the right 
BPS energy bound, expected of fundamental strings orthogonal toaD(p-l) brane. 

7.1 Smooth Domain Walls Threaded by String Fluid 

When V is included in the analysis, we may consider a different kind of static so- 
lutions, involving domain walls. A singular solution of domain wall type has been 
known M, but recently it was found that a smooth domain wall solution is possi- 
ble when we modify the asymptotic boundary condition [[H], |32|, B3J. Here we will 



reproduce these solutions within our formalism and study its properties. 

For this, we will consider static configurations with Pi = and ttt = 0, and 
assume a uniform electric flux, ir x . We will achieve Pi = by setting Fij = 0. With 
this, another allowed quantity is 

Vt = n x d x T. (7.1) 
In such configurations, a useful identity is 



h = vw + v(m(i + (d x m. (7.2) 



Then the momentum conservation (6.7) forces 



which may be rewritten as 

C 2 (ir 2 x + V(Tf) = 1 + (d x T) 2 (7.4) 

for some constant C. Inspection of (|6.1|) shows that this constant is related to the 
asymptotic value of E x — > e as (Ctt x ) 2 = e 2 .^ Thus we only need to solve 

% = ±vW + W!R- (7-5) 

There are three generic cases to consider 

TTFor the present solutions E x itself is a constant, by the way. 



- 26 - 



• If e 2 is larger than or equal to 1, this gives monotonic solution T(x), and 
represents a single D(p — 1) (or anti-D(p — 1)) brane at T = 0. 

• If e 2 is less than 1 but larger than 1/(1 + V(0) 2 /7r 2 ), this gives an oscillatory 
T{x) with an infinite array of zeros and represents an alternating array of 
D(p — 1) and anti-D(p — 1) branes. 

• If e 2 is less than 1/(1 + V^O) 2 /^ 2 ), no solution exists. 

In all cases where solutions exist, the configuration is that of smooth domain wall 
solutions, generically threaded by uniform distribution of flux lines transverse to the 
domain wall. 

7.2 A Minimal Solution and BPS Energy 

Energy density of the above solutions are 



H = yM+VW+Wm = + V 2 ) (7-6) 

7T~ 



which approaches en x asymptotically. We are interested in solution with pure string 
fluid far away from domain wall, and this forces us to consider e = 1, in particular. 
This choice is also energetically favored since it represents the lowest energy solution 
available, given the conserved flux and the single domain wall. As we will see shortly, 
this case may be regarded as a BPS-saturated solution. 

As would be expected from — 1, the equation of motion simplifies to 

£ = ±332 (7.7) 
ax ti x 

whose solution is such that T does approach vacuum at x — ±oo. The energy density 
also simplifies further for e = 1 as**, 

H=\n x \ + ^r = n x + V\d x T\ (7.8) 
Fx | 

The first piece is the BPS energy density associated with the fundamental string 
charge. The second piece must be associated with the domain wall itself. In fact, 
integrating over x, this is precisely the tension of the D(p — 1) brane 

Jvd x Tdx = J VdT = r p _ 1 (7.9) 

The energy of the solution is then precisely sum of two terms, one from tension 
of the D(p — 1) brane and the other from transverse fundamental string. This gives 



"We are indebted to Chanju Kim on this point. 
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exactly the BPS energy expected of fundamental string ending on or passing through 
a D-brane. 

Thus at least when we have uniform string fluid threading a D(p- l)-brane, 
the string fluid can mimic, quantitatively, behavior of fundamental string ending on 
D-branes. This may be compared with findings of Ref. [f£|], where it is asserted that 
string fluid cannot end on D(p — 1) brane. The latter is based on a computation 
where an infinitely thin D(p — 1) is taken to be a background and neglects a possible 
backreaction of the domain wall solution to the presence of string fluid. While the 
assertion might stand when transverse string flux are well-isolated, this example 
cautions us to be careful about backreaction of a domain wall. 



8. Summary 

We have reviewed dynamics of the tachyon coupled to a gauge field. When the 
tachyon condensation has progressed far so that V — > 0, the dynamics is that of 
two fluids, string fluid and tachyon matter, as anticipated, but the fluid equation of 
motion must be augmented by a set of integrability condition, which is necessary if 
we wish to recover and T from the fluid variables. We have isolated a large family 
of static solutions with huge degeneracy and tested their stability. 

We further extended the formalism to the case of V ^ somewhere, and showed 
that during initial stage of condensation, electric flux lines tends to be attracted 
toward T = 0. Since T = may survive the condensation process only if there is a 
domain wall, this initial configuration will lead to a D(p — 1) brane with fundamental 
string flux spread on it. Interaction between a domain wall and transverse string 
fluid is more drastic, and the latter thickens the former. We also wrote down such 
smooth domain wall solution, using the Hamiltonian formulation, and discovered that 
minimal solution of this type saturates a BPS bound of strings ending on D-branes. 

One aspect of the V — limit, worth emphasizing, concerns the coupling between 
the two fluid components. Despite a rather tight coupling between the string fluid and 
the tachyon matter, as evidenced by the form of the Hamiltonian, static distribution 
of the two fluid components seems pretty much independent of each other. The only 
real constraint is that static distribution of tachyon matter must be uniform along 
the string fluid direction. One might have expected that two fluid components come 
with the same (transverse) distribution, but this is not the case at all. Instead, the 
tight coupling between the tachyon matter and the string fluid affects the subsequent 
dynamics in an unexpected way. As is clear from the fluid equation, distribution of 
string fluid and distribution of tachyon matter evolves in the complete absence of 
pressure (except for the tension along the string fluid) and is essentially free. That is 
they simply responds to velocity field of fluid, which is in turn affected by gradient 
of density of flux direction along the flux lines. Regardless of the details of the latter, 
however, two fluid shares one and the same velocity field v. This implies that once the 
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relative distribution of the two fluid component is set, the subsequent perturbation 
moves string fluid and tachyon together in such a manner that keeps this ratio fixed 
when followed by co-moving observers. 

The classical system of tachyon and gauge field turned out to have many intrigu- 
ing surprises; fluid-like behavior of the perfectly sensible field theory, the drastic 
reduction of perturbative degrees of freedom, Carollian collapse of effective light- 
cone, and finally many intriguing classical solutions, static or dynamic, singular or 
nonsingular. One question that remains murky is whether this classical field theory 
admits a sensible and practical quantum treatment of its own, regardless of its tie to 
stringy context, and whether interesting physics arise from a quantization. We feel 
that this avenue of research deserves more study in future. 
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A. S-Dual Formulation of V = Limit 



Ref. [£0J discovered a nonvanishing dual Lagrangian for pure gauge system, whose 
equation of motions and Bianchi identity are respectively Bianchi identity and equa- 
tion motion of the original Born-Infeld system. It is rather clear that dualization 
procedure works even when scalar fields are present, once we think of the latter as 
components of the gauge field along some fictitious dimensions. 

The prescription for the dualization process is to introduce a new gauge field C 
and its field strength, G = dC such that 

/C = 7r m dt A dx m + — K mn dx m A dx n = *G. (A.l) 

Momentum ir m is canonical conjugate momenta of A m = [A±, • • • , A p , T). K mn are 
defined as 

K_-g-. (A,) 

u± mn 

Finally, the Hodge-dual operator * here is defined with respect to the "volume form" 

dt A dx 1 A dx 2 A ■ • • A dx p A dx T . (A. 3) 

With this, the Legendre transformation of the Hamiltonian gives a dual Lagrangian 
as a function of the dual field strength G, 

C! = V-/C72 = y/(P/2. (A.4) 
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One crucial point here is that the Legendre transformation that leads to £' is not 
reversible. Thus, to obtain truly dual Lagrangian, one must introduce constraints on 
the dual side. 



As explained in Ref. [20], this does not complete the dualization process due 
to the fact that the original Hamiltonian is quite degenerate. That is, part of F mn 
which has no inner product with 7r m never enters the Hamiltonian. Because of this, 
inverse Laplace transform of £' does not give us back TC. This problem is curable by 
adding constraints, 

K A K = 0. (A.5) 
Thus the correct dual Lagrangian may be written as 



\J G 2 /2 + (A, *G A *G) (A.6) 
with a Lagrange multiplier field A. This is where the pure gauge dynamics maps to 



that of Ref. [45 , which attempted to write a field theory that produces string-like 



degrees of freedom classically. 

B. Canonical Field Equation with V = 

The tachyon effective action with gauge field is given by 

S= f d p+1 x £, (B.l) 



where 

£ = -V(T)\/=X, 

= 9jw + Ffa, + d^TdvT, 
X = detX^. (B.2) 

Then the determinant X can be expressed as 

X = Aoo detpfy) - XuD^Xjo, (B.3) 

where the matrix D is transpose of the cofactor for matrix (X)ij and the components 
of matrix {X)^ v are written by 

X 00 = —N 2 + h l3 VV + T 2 , 
Xni = Et , 



Xio = — E, 



— hij + Fij + diTdjT 

with 



Ef = F 0i ±Td i T±L"h ji . (B.4) 
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Let us denote the Hamiltonian using canonical variables to describe the dynamics 
of the system by Hamiltonian equations. Then conjugate momenta are given by 

. <)£. V EfD ji + D ij Ej 

7T = — - = — — , (B.5) 

8C V (■ EfDi%T - diT&i E- 
n T = — = -p= I TdetXij - 3 - 



df 



(B.6) 



where n l (tt t ) is conjugate to Aj (T), and n % satisfies the Gauss constraint = 0. 
The Hamiltonian is obtained by the following Legendre transformation 

H = n'Ei + ix T f-C + 7r l diA . (B.7) 

From now on, let us use matrix notation for simplicity. Then the quantities which 
we have defined are denoted by in matrix forms in temporal gauge A Q = 0, 

X = X 00 det Xij + E + DE~, (B.8) 

V E+D + DE- , , 

* = 7=7 2 • < M > 

^.^(fd^- ggg^ggg: ), ( B,o, 

= V } E + DE + EDE- f2 _ E+D8TT - TdTDE' 

V^X ^2 3 

V , s , E + DhL- LhDE-\ 

= -j= [{N 2 - LhL) det ^ + \ , (B.ll) 

where all matrix indices are k = 1, • • • ,p. Let us define a matrix for convenience 
in calculations 

X = h + F + 8T8T, (B.12) 

which has properties 

DX = XD = det XijI, (B.13) 
where J is p x p unit matrix. Then we obtain the following relations 

vrXvr = tt/itt + (ndT) 2 , 

V f^^, „ E+-E-, E+Dh-hDE-\ 

Fix + dTn T = —= TdT det X^ det + . 

y—X \ 2 2 J 

(B.14) 

Using the relations in Eq. ( |B.14|) , we find 

V_ 
y/^X 



nhn + 4 + (Ftt + dTTX T )h- l (Fit + dTn T ) + (ndT) 2 + V 2 det X {j = -j=N det X {j , 



(B.15) 
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where the matrix h 1 has components (h = . Now we can rewrite the Hamil- 
tonian representation using canonical variables as 



n 



V 



N 2 det X i:j - ttFL + ix T dL 



N^/nlm + Trf, + (Fvr + dTic^h-^Fir + dTn T ) + (vrdT) 2 + 4 + 1/2 det x .. 
—ttFL + tttOTL, (B.16) 



where we used the relation 



E + DLh - LhDE- = 2 \LhL det Xy - ^—^-nFL + ^-^-7r T <9TL^ . (B.17) 



Then the Hamiltonian equations in V = limit are given by 
f 



7Tt 



dU 


N 


dn T 


W 


on 




dT 




m 


N 


drc 1 


Vy 



(tt t (1 + OTh^dT) + dTh^Fn) + dTL 
A (7r l (7r9T) + TT T h ij {F jk 7i k + ^Tttt)) 



(B.18) 

+ 9 i (7r T i i ), (B.19) 
(FL)i, (B.20) 



7T 



^ « ( N 



Y 



y k - 7c j h lk )(Fn + dT7i T ) k ^j - djtfU 



where Y = irhn + vr 2 + (Fvr + dTir^h-^Fir + <9Tvr T ) + (vr<9T) 2 . 



(B.21) 



C. Second Order Perturbation Equations 

General perturbation equations for Eqs. (|5.1|), ( p.3[ ) are given by, 



2n 



O.r 



1 71—1 



5+ 6 » = 



i=l 

71-1 



2no x y 2n 0x 



i=l 



7i-l 



n-l 

2n(o).6H = _^6».^-). 



eft* ■a^ n - i \ 



(C.l) 
(C.2) 
(C.3) 
(C.4) 



i=l 
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Now let us consider the second order perturbation equations. The y-components 
of the second order perturbation equations in Eqs. (|C.1|) , (|C.2|) are written by 



2n, 



-(& (1) -<9)a« 



Ox 



2n 0x {2n 0x ) z J 2n 0x 2n 0x 



(C.5) 



2n 



Ox 



2n 0x {2n 0x ) z J 2n 0x 2n 0x 



where we used the results of the first order perturbation ( J5.8[ ) and (5.10). Then we 
obtain the solutions of these inhomogeneous first order differential equations, 



,( 2 ) 



a v-v = h(x + ,y) 



1 



n: 



G + ^Lt F + F) + 



2n 0x 2n 0x 



(2n 0a! 



n X jt\ j(l,0) ' 



2n, 



-Ff 



(o,i) 



Ox 



bP = h( X -,y) - ( J-G - ^ (F + F) + ^jf) /('■») - -J-F/M, 
\zn 0x lriQ X v / \lnQ X ) J zn 0x 



(C.7) 



(C.8) 

where h and ft, are arbitrary functions and we define, 

/x + p X ~ 

F(x + ,y) = dwF(w,y), F(x-,y)= dw F(w,y). (C.9) 



Using the Eqs. (|5.8|), ( |5.10| ), ( |C.7| ), we can obtain the second order perturbation 
equations for the x-components of the Eqs. ( |C.1| ),( P72] ) as follows: 



2n 0x V y 



2n 0x (2n te ) 2 " 



2n 0x (2n 0x y 



n o x 
2n 0x 



tf)g 



,(1,0) 



+ 
+ 



n 0x £f _|_ ( n 0z) 2 



n; 



/ \2 
Ox J 



(2n te ) 2 



(2n 0:r ) 



(2n 0a: ) 



t F + F / 



?(1,0) 



r?. 



O.r 



n; 



Ox) 2 



2 K(d 2 n te ) - «) 2 ) - « ^ 



1 (0,1) , n 0x £(0,1) 



(2n 0x y 9 (2n 0x ) 3 ~ 2n 0x 1 ' ' (2ra x' 12 



''OxV 



^6l 2 ) 



(2n x) 3 
4 2) < + " W 



(2n 0x y 



tf)f, 



(CIO) 



2n, 



Ox 
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2n 0x {2n 0x ) 2 \ 2n 0x {2n 0x ) 2 2n 0x J 

, ( n 'ox q , KJ jr I KJ + ( F + p\\ f(ifi) 
+ V(2n ,) 2 + (2n ,) 3 + (2n 0x f V ) J 1 

( n 0x _ ( n 0x) 2 p j_s(0,l) , n 0x F (0,l) 

\(2n 0x y y {2n Qx f 2n 0x y (2n 0x ) 2 



2 {n 0x (d 2 n 0x ) - (n' 0x ) 2 ) 



In: 



I \2 

Ox 



(2n 0x ) 3 ' (2n 0x ) 

Using these results, we find 

a a „(2) _ ( n Qx _ ( n 'ox) 2 p K) 2 ff \ 1(1,0) 

+ ^ - \{2n ax y 9 {2n,xf (2n 0x f J J 1 



F +77^tf\f- (Cll) 



n 0x 



^ //(0 ' 1,+ ^ (/ ;^ (1 ' 0, -i /9 7 

+ \ ———g + — -F - — ^-i/ j g (2 ' 0) - (/ + /)/ 

V 2n 0a; (2n x) 2 2n 0x / (2n 0:E ) 3 

(C.12) 



I n to f f (0,1) _ n o* / f f^(l,0) ^f^ 1 ' 1 ) 

+ (2n ,) 2// (2n ,) 2lJ JJy 2nJ 9 

+ (~9 + j^t-J + ^Ltf) - p^-Jf + f)f, 



2n 0x (2n 0x ) 2 2n 0x J (2n 0x ) 3 

(C.13) 

(2n 0a; ) 2 2n 0a: 

1 + f) /( 2 -°), (C.14) 



2n 0x 2n 0x (2n 0x ) 
y (2n 0a; ) 2 2n 0x 



2n x 2n 0:r (2n O a0 
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